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Abstract 

We consider a general N-degree-of-freedom nonlinear system which is chaotic 
and dissipative and show that the nature of chaotic diffusion is reflected in the 
correlation of fluctuation of linear stability matrix for the equation of motion 
of the dynamical system whose phase space variables behave as stochastic 
variables in the chaotic regime. Based on a Fokker-Planck description of the 
system in the associated tangent space and an information entropy balance 
equation a relationship between chaotic diffusion and the thermodynamically- 
inspired quantities like entropy production and entropy flux is established. 
The theoretical propositions have been verified by numerical experiments. 
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I. INTRODUCTION 



Several authors have enquired recently about the relationship between phase space dy- 
namics of a dynamical system and thermodynamics The question acquires a special 
relevance for the dissipative system when the phase space volume contracts by virtue of 
possessing the attractors and also when the system is nonlinear and comprises a few-degrees- 
of-freedom. Thus even when these systems are not truly statistical in the thermodynamic 
sense, it is possible that chaotic diffusion due to intrinsic deterministic chaos or stochastic- 
ity plays a significant role in the dynamics. It is therefore worthwhile to enquire about the 
relationship between chaotic diffusion in a dynamical system and the thermodynamically- 
inspired quantities like entropy production and entropy flux. Our purpose in this paper is to 
address this specific issue. 

In what follows we shall be concerned with the nonlinear dynamical systems which are 
chaotic and dissipative. We do not consider any stochasticity due to thermal environment 
or external nonthermal noise. The "deterministic stochasticity" (i.e, chaos) has a purely 
dynamical basis and its emergence in nonlinear dynamical system is essentially due to loss 
of correlation of initially nearby trajectories. This is reflected in the linear stability matrix 
or Jacobian of the system When chaos has fully set in, the time dependence of this 
matrix can be described as a stochastic process, since the phase space variables behave as 



stochastic variables [111. It has been shown that this fluctuation is amenable to a theoretical 



description in terms of the theory of multiplicative noise |T^. Based on this consideration 
a number of important results of nonequilibrium statistical mechanics, like Kubo relations, 
fiuctuation-decoherence relation, fluctuation-dissipation relation and exponential divergence 
of quantum fluctuations have been realized in chaotic dynamics of a few-degree-of-freedom 
system |T^-|18[. In the present paper we make use of this stochastic description of chaotic 
dynamics to formulate a Fokker-Planck equation of probability density function for the 
relevant dynamical variables, the "stochasticity" (i.e, chaoticity) being incorporated through 
the fluctuations of the time-dependent linear stability matrix. Once the drift and chaotic 
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diffusion terms are appropriately identified the thermodynamic-like quantities can be derived 
with the help of the suitable information entropy balance equation. 

The outlay of the paper is as follows: In Sec. II we introduce a Fokker-Planck description 
of the dynamical system and identify the chaotic drift and diffusion terms. This is followed 
by setting up of an information entropy balance equation in Sec. III. We then look for the 
entropy flux and entropy production-like terms in the steady state. The shift of stationary 
state due to additional external forcing and the associated change in entropy production is 
considered in Sec. IV. We illustrate the theory in detail with the help of an example in Sec. 
V. The paper is concluded in Sec. VI. 



II. A FOKKER-PLANCK EQUATION FOR DISSIPATIVE CHAOTIC 

DYNAMICS 

We are concerned here with a general N-degree-of-freedom system whose Hamiltonian is 
given by 

N 2 

^ = E^ + nte}.^) . i^l---N (1) 

where {qiiPi) are the co-ordinate and momentum of the i-th degree-of-freedom, respectively, 
which satisfy the generic form of equations 

dH 

OPi 

dqi 

We now make the Hamiltonian system dissipative by introducing —■jpi on the right hand 
side of the second of Eqs.(2). For simphcity we assume 7 to be the same for all the N degrees 
of freedom. By invoking the symplectic structure of the Hamiltonian dynamics as 

Qi ior i — 1- ■ ■ N , 
Pi-N iori^N + l,---2N , 

V 

and defining I as 



E 
-E --fE 

where E is an iV x unit matrix, and is an x null matrix, the equation of motion 
for the dissipative system can be written as 



. dH 

j=l ^^j 



(3) 



We now consider two nearby trajectories, Zi, Zi and Zi + Xj, Zi + Xi at the same time t in 
2N dimensional phase space. The time evolution of separation of these trajectories is then 
determined by 

2N 

Xi = Y^ jij(t)Xj 

in the tangent space or separation co-ordinate space Xi, where 



(4) 



Jij — X/ 



ik 



dzkdz^ 



Therefore the 2N x 2N linear stability matrix J assumes the following form 



(5) 



J 



(6) 



E 

M{t) -7^ 

where M is an x matrix. Note that the time dependence of stability matrix J_{t) is 
due to the second derivative -^§^ which is determined by the equation of motion (3). The 
procedure for calculation of Xj and the related quantities like Lyapunov exponents is to solve 
the trajectory equation (3) simultaneously with Eq.(4). Thus when the dissipative system 
described by (3) is chaotic, J(t) becomes ( "deterministically" ) stochastic phase space due to 
the fact that Zi-s behave as stochastic phase space variables and the equation of motion @ 
in the tangent space can be interpreted as a stochastic equation [p!3H T8|. 

In the next step we shall be concerned with a stochastic description of J_{t) or M{t). For 
convenience we split up M into two parts as 



M = Mo + Mi{t) 



(7) 



where Mq is independent of variables {zi} and therefore behaves a constant part and Mi is 
determined by the variables {zi} for i = 1 ■ ■ ■ 2N. Mi refers to the fluctuating part. We now 
rewrite the equation of motion (4) in the tangent space as 

X = JX 

= L{{X,},{z,}) (8) 

where X and L are the vectors with 2N components. Corresponding to L in (|]) can be 
split up again to yield 

X = L%X) + L\X,{z,{t)}) 
or X, = L,\{X,}) + L]{{X,},{z,}) , t = l---2N. (9) 

Eq.(^) indicates that Eq.(|D is linear in {Xi}. Eqs.(^), (||) and (||) express the fact the 
first N components of are zero and the last components of are the functions of 
{Xi} ioT i = 1 ■ ■ ■ N only. The fluctuation in Lj is caused by the chaotic variables {zi}-s. 
By deflning Vx as differentiation with respect to components of X, i. e. , {Xi} (explicitly 
Xi = AqifoTi = l---N and Xi = Api for i = N + 1 ■ ■ ■ 2N) and since L} = for i = 1 ■ ■ ■ N 
and L} = Li (Xi ■ ■ ■ Xjv) for z = iV + 1 ■ ■ ■ 2N we have V x ■ = Etii^ ' ^) + S'fiv+i ^ ' 
L}{Xi ■ ■ ■ Xn) = 0. This allows us to write the following relation (which will be used later 
on), 

Vx ■ L'mX^}) = ■ Vx<P{{X,}) (10) 

where 0({Xj}) is any function of {Xi}. 

Note that Eq.(^ by virtue of (|]) is a linear differential equation with multiplicative 
"noise" due to {zi} determined by equation of motion (|). This is the starting point of our 
further analysis. 

Eq.(|^) determines a stochastic process with some given initial conditions {Xi{0)}. We 
now consider the motion of a representative point X in 2N dimensional tangent space 
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(Xi ■ ■ ■ X2n) as governed by Eq.(|]). The equation of continuity, which expresses the con- 
servation of points determines the variation of density function t) in time as given 
by 

d(t>{X,t) 



dt 

Expressing Aq and as 



-V^-L{t)<f){X,t) . (11) 



A, = -Vx ■ 



10 — - V X 

and = -Vx ■ (12) 

we may rewrite the equation of continuity as 

^^i^ = [Ao + aA,itMX,t) . (13) 

It is easy to recognize that while Ao denotes the constant part Ai contains the muhiphcative 
fluctuations through the phase space variables of the dynamical system {zi{t)}. a is a pa- 
rameter introduced from outside to keep track of the order of fluctuations in the calculations. 
At the end we put a = 1. 

One of the main results for the linear equations of the form (|T3|) with multiplicative noise 



may now be in order ||T2|. The average equation of {(f)) obeys [ P{X,t) 



P = !^Ao + a{Ai) dT{{Ai{t)exp{TAo)A,{t-T)))exp{-TAo)^P{X,t) . (14) 

The above result is based on second order cumulant expansion and is valid when fluctu- 
ations are small but rapid and the correlation time Tc is short but finite or more precisely 

{{A,{t)A,{t'))) = for \t-t'\ > r,. (15) 

We have, in general, (A-^) ^ 0. Here ((■■■)) implies ((OCj)) = idCj) - {QiCj)- 

The Eq.(|T^) is exact in that limit Tc 0. Making use of relation (|T2|) in ( [TT| ) we obtain 



dP 
'dt 



|-Vx ■ L° - a{Vx ■ L') + a' rfr((Vx ■ L\t) exp(-rVx ■ L°) 

Vx ■ L\t - r))) exp(rVx ■ L^)} P . (16) 



The above equation can be transformed into the following Fokker-Planck equation (a = 1) 
for probabihty density function P{X,t), (the details are given in the Appendix A); 



dPjX, t) 
dt 



dXidXj 



(17) 



where, 



F = L° + (L^) + Q 



(18) 



and Q is a 2A^-dimensional vector whose components are defined by 



{{R'.))dTdeti{T)det2{T) 



(19) 



Here the determinants deti, det2 and i?' are given by 



detiir) 



dX-^ 



dX 



and det2{T) 
d 



dX 



dX- 



and R'^=J2Ll{X,t)^Y.LliX-\t-T) 



« k 



dX^^ 



(20) 



It is easy to recognize F as an evolution operator. Because of the dissipative perturbation 
we note that div F < 0. 

The diffusion coefficient Vij in Eq. ([T7|) is defined as 



T.{{Ll{X,t)Ll{X-\t 



dX. 

] 

dXl 



-))dr 



(21) 



k "'^'■k 

We have followed closely van Kampen's approach [|12[ to generalized Fokker-Planck equation 
(p!7|). Before concluding this section several critical remarks regarding this derivation need 
attention: 

First, the process Mi(t) determined by {zi} is obtained exactly by solving equations of 
motion (^ for the chaotic motion of the system. It is therefore necessary to emphasize that 
we have not assumed any special property of noise, such as. Mi (t) is Gaussian or 5-correlated. 
We reiterate Van Kampen's emphasis in this approach. 

Second, the only assumption made about the noise is that its correlation time Tc is short 
but finite compared to the coarse-grained timescale over which the average quantities evolve. 



Or in other words the velocity changes should be small, smooth and uncorrected after short 
times. This assumption, however, puts a restriction on the applicability of the present theory 
to certain class of systems, e. g. , systems subjected to 'hard' collisions such as billiards 
and molecular systems also in certain non-dissipative Hamiltonian systems like standard 
map for which the usual assumptions about the rapid decay of correlations and fluctuations 
are not valid and entropy production does not occur. The special reference may be made 
in this connection to the work of Zaslavsky and collaborators ||19|] to demonstrate that 
in many real systems the decay of correlation exhibits power law dependence, distributions 
admit of infinite moments and the fluctuations become long-lasting. Since the mathematical 
difficulties of dealing with finite arbitrary correlation time of noise in a chaotic system is 
quite formidable we confine ourselves in the present discussion to chaotic systems with short 
but finite noise correlation time. 

Third, we take care of fiuctuations upto second order which implies that the deterministic 
noise is not too strong. 

Eq.([T7|) is the required Fokker-Planck equation in the tangent space {Xj}. However the 
important point is to note that the drift and diffusion terms are determined by the phase 
space {zi} properties of the chaotic system and directly depend on the correlation function 
of the fluctuations of the second derivatives of the Hamiltonian (15|). 



III. INFORMATION ENTROPY BALANCE: ENTROPY PRODUCTION 

We shall now consider the well-known relation between probability density function 
P{X, t) and information entropy S as given by 

S=- J dXP{X, t) In P{X, t) . (22) 

Note that in the above deflnition of entropy we use P{X,t), the probability distribution 
function in the tangent space since one is concerned here with the expansion of the phase 
space in terms of a tangent space and dialation coefficients of the dynamical system for which 



8 



the expanding and contracting manifolds can be defined. On tlie otlier hand it is worthwhile 
to recall the dynamic entropy of a dynamical system (Kolmogorov entropy) defined in terms 
of the properties of evolution in the tangent space. A remark on the connection between 



entropy and expansion by Sinai ||20| is noteworthy in this context ; " It already seems clear 
that positiveness of the entropy and presence of mixing is related to extreme instability of 
the motion of the system : trajectories emanating from the nearby points must, generally 
speaking, diverge with exponential velocity. Thus entropy is characterized here by the speed 
of approach of the asymptotic trajectories" which is formalized by defining the expansion 
coefficient as logarithm of the relative increase under the flow of a volume element in the 
expanding manifold. Our definition of information entropy (22) makes use of the tangent 
space description of the systems in terms of logarithm of the probability of expansion in the 
tangent space, keeping in mind that — In P is "a measure of unexpectedness of an event (the 
amount of information) and the information entropy is a mean value of this unexpectedness 
for the entire system" |21|]. The definition (22) is therefore different from Kolmogorov 



entropy. We emphasize that even in absence of any direct formal connection between P{X, t) 
and the phase space distribution function it is possible to use the distribution function 
P{X,t) defined in the tangent space to have an explicit expression for entropy production- 
like quantity as a function of the properties of phase space variables {zi} of the dynamical 
system (i. e. , in terms of drift and diffusion coefficients of the Fokker-Planck equation). 
The above definition of an information entropy like quantity allows us to have an evolu- 



tion equation for entropy. To this end we observe from Eqs.(|T^ and (0) that p^j23 



dS_ 
'dt 



InP . (23) 



Note that the probability density function P(X,t) is defined in tangent space {Xi\. V 
and F as expressed in Eqs. (pi]) and (|18|), respectively, are determined by the correlation 
functions of fiuctuations of second derivative of the Hamiltonian of the system. Eq.(p3D 
therefore suggests that the entropy pro duct ion- like term originating from Eq. ( ^5] ) is likely to 
bear the signature of the chaotic dynamics. The relation is direct and general as it is evident 
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from the following equation [obtained after partial integration of Eq.(^) with the natural 
boundary condition on P{X, t) that it vanishes as |X| oo and assuming the X-dependence 
of Vij to be weak (as a first approximation)]; 

ds r ^ w-D [ ^ dP dP 

The first term in ( ^4]) has no definite sign while the second term is positive definite 

as 



f ^/«Pv..f+i:i:D.,/i|£|£« . (24) 



because of positive definiteness of Pjj. Therefore the second one can be identified p2 
the entropy production 

S.o.^EEA./^llll^ (25) 

in the steady state. The subscript s of Pg refers to steady state. It is evident from Eq.(plD 
that 

Sfi^, = jdX P,(X) Vx-F = VlTF 

Sprod = —Sflux ■ (26) 

Note that since the chaotic system is dissipative Vx-P is negative [See Eq.(|^)]. 

It is thus evident that the relations (^), (^) illustrate the dynamical origin of entropy 
production like quantity in a chaotic dissipative system. The dynamical signature is mani- 
fested through the drift term P and the chaotic diffusion terms in Vij. It must be emphasized 
that the notion of diffusion has nothing to do with any external reservoir. Rather it pertains 
to intrinsic diffusion in phase space of the chaotic system itself. 



IV. THE CHAOTIC SYSTEM DRIVEN BY AN EXTERNAL FORCE 

We shall now examine the entropy production when the dissipative chaotic system is 
thrown away from the steady state due to an additional weak applied force. To this end we 
consider the drift Pi due to external force so that the total drift P has now two contributions 

P{X) = Po{X) + hPi{X) . (27) 
10 



When h = 0, P = Ps- The deviation of P from Pg in presence of nonzero small h can be 



explicitly taken into account once we make use of the identity for the diffusion term [52 



d^P 



d 



P 



(91nR 



+ 



d 
dX, 



P. 



d P 

djqp's 



(28) 



When P = P^ the second term in 



vanishes. In presence of additional forcing the 



Eq.(|17D becomes, 

dP _ 
'dt~~ 

where if) is defined as 



-Vx.^P - hVx.FiP + E E 



d 
dX, 



P. 



d P^ 

djqp:^ 



(29) 



dlnPs 
dXi 



(30) 



Here we have assumed for simplicity that Vij is not affected by the additional forcing. The 
leading order influence is taken into account by the additional drift term in Eg. (p9|) . 

Under steady state condition (P = P^) and h = 0, the second and the third terms in 
(53) vanish yielding 



Vx.i^Ps = o . 



(31) 



It is immediately apparent that ipPs refers to a current where J' = ipPg- The steady state 
condition therefore reduces to an equilibrium condition {J = 0) if 



^ = . 



(32) 



(In Sec. V we shall consider an explicit example to show that = 0). This suggests a formal 
relation between Foi and Vij as 

d\nP, 

3 



(33) 



where Ps may now be referred to as the equilibrium density function in separation co- 
ordinate space. Fq contains dissipation constant 7 and the diffusion matrix T>ij is a function 
of correlation function of fluctuations of the second derivative of the Hamiltonian. 



11 



To consider the information entropy balance equation in presence of external forcing we 
first differentiate Eq.(E2[) with respect to time and use Eq.(E9]) to obtain 



^5" d , ^ 
— = -- / dXPlnP. 
at at 



^ 3 



3 s , 



(34) 



It is apparent that as P deviates from P,, differs from unity and the entire second 
integral within the parenthesis is of the second order. Note that InP^ in the first 



integral in Eq. ([3^) is a constant of motion and the integral denotes its average. The first 
term vanishes because it is of higher order as it involves P^ and others. (Moreover since in 
the discussion that follows we consider the steady state, this term does not contribute to 
the subsequent calculations). To compute the contribution AS" to entropy balance due to 
the external forcing only we perform integration of the second, third and fourth terms by 
parts. We thus obtain, 

^^,./..wv.....,./.x(i:...^?|^).P 

+ (35) 

where we have put h5P = P — Pg. In the new steady state (in presence oi h ^ 0), the 
entropy production and the flux like terms balance each other as follows: 

^Sprod = —^Sfiux (36) 



with 



^S^rod = y Vi. [ dXP ( In — I ( In — 

'J \dX, Ps) \dXj Ps, 



(37) 



and 



^S^iux = h'J dX6PVx.Fi + h'JdX ^1^^^) 



(38) 



In the following section we shall work out a specific example to provide explicit expres- 
sions for the entropy production and some related quantities due to external forcing. 
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V. APPLICATIONS 



A. Entropy production in the steady state 

To illustrate the theory developed above, we now choose a driven double-well oscillator 
system with Hamiltonian 



H = + aq^ — bqf + eqi cos ilt 



(39) 



where pi and qi are the momentum and position variables of the system, a and b are 
the constants characterizing the potential, e includes the effect of coupling constant and the 
driving strength of the external field with frequency Q. This model (|39|) has been extensively 



used in recent years for the study of chaotic dynamics [13-15,2^. 

The dissipative equations of motion for the tangent space variables Xi and X2 corre- 
sponding to qi and pi (Eq.8) read as follows: 



d 
di 



where J as expressed in our earlier notation 

















= J 




■1 




X2 




X2 






= X2 J 



zi = qi Z2= pi 



is given by 



J 



^ £ ^ 



-E --fE 



d'^H a' 



dzidzi dzidz2 
\ dz2dzi dz2dz2 / 



Therefore J reduces to 



1 \ 

C(i) + 2t -7^ 



where ({t) = —12azf. Thus we have 



(40) 



Mo = 26 Ml = C(t) 
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Eq.(^OD is thus rewritten as 



d 
It 



(41) 



with 



X2 

2bXi - 7X2 



and 



where L° and are the constant and the fluctuating parts, respectively. The fluctuations in 
L^, i.e., in ((t) is due to stochasticity of the following chaotic dissipative dynamical equations 
of motion; 



Zl = Z2 

Z2 = —azf + 2bzi — e cos Qt — 72:2 . 
Now for the constant part and the fluctuating part we write 



(42) 



L"' = X2 L"^ = 2bXi - 7X2 
= L^^ = C(t)Xi . 

We may then apply the result of Eq. ( |A5D . 

The mapping X X* is found by solving the 'unperturbed' equations 

Xi = X2 

X2 = G2 — 7X2 

Comparison with Eq. shows that 6*2 (= 2bXi) is free from X2. 

As a short time approximation we consider the variation of Xi and X2 during ; 

Xr = -rX2 + Xi = Gi(Xi,X2) 

= -TG2 + e^^X2 = 6-2 (Xi, X2) . 

So the 5f-matrix of Eq. ( |A15|) becomes 



(43) 
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^ 1 re"^^ ^ 



The vector R from Eq.( |A15|) can then be identified as 



R 



U{t-r){Xi-rX2)gi2^ 



at-r){X,-rX,)g22 

r \ 



I 



at - r)Xire-^- 
C(t-r)(Xi-rX2)e-^- 

Similarly the vector R' is given by 



neglecting the terms of 0(7 



R' 







-at - r)re-^-XiC(t) 



From Eqs.(43) and (44) we have 



deti{T)det2{T) ~ 1 



Then the vector Q can be written as 

/ 

Q = 







^Xi/o-((C(t)C(t-r)))re-^^rfr ] 
Now the diffusion matrix T> can be constructed as 



V 



S ^ 



V 



P2I T^22 



where 



V2i = xm / mt)at-r)))re-^^dT 

Jo 



and 



= XtiO) / ((C(t)C(t - r)))e-^^rfr - X^MO) / ((CWC(t 
JO Jo 
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It is important to mention that the assumption of weak X-dependence of the diffusion 
coefficient (by freezing its time dependence) is permitted as a first approximation within the 
perview of the present second order theory for which the strength of noise is not too large. 
We also emphasize that for actual theoretical estimate of the entropy production in terms of 
the formulae (26) or (58) explicit evaluation of the diffusion coefficients is not required (see 
the next section). A straightforward calculation of drift is sufficient for the purpose. This 
point will be clarified in greater detail in the next Sec.VB. 

Then the Fokker-Planck equation (17) for the dissipative driven double- well oscillator 
assumes the following form: 

dt dXi ^dX2 'dX2^ ' 3X28X1 ^^dXl ^ ' 

where 

a;^ = 26 + c + C2 

with 



00 



C2= / {mat-r)))Te-^'dT 
Jo 

c = (C) . (51) 



The similarity of the equation (50) to generahzed Kramers' equation can not be overlooked. 
This suggests a clear interplay of chaotic diffusive motion and dissipation in the dynamics. 
We now let 

U = asXi + X2 (52) 

where is a constant to be determined. 

Then under steady state condition Eq.(50) reduces to the following form : 

A(A.t/)A + P.^ = , (53) 

where 
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'Ds = T>22 + asT>2i , 

and 

XM = -asX2 - up-Xx + 7X2 . 

Here Ag is again a constant to be determined. Putting (52) in A^C/ as given above and 
comparing the coefficients of X\ and X2 we obtain 

AgOs = —up' and A^ = — + 7 . 

The physically allowed solutions for ag and As are as follows ; 

7 - V7^ + 4u;^ 7 + Vt^ + 4a;^ /..x 
Os = and As = . (54) 

The stationary solution of (53) Pg is given by 

P^^Ne"^ . (55) 



Here N is the normalization constant. By virtue of (55) t/j corresponding to Eq.(30) is 
therefore 

V^ = Asf/-Ps^^ = . (56) 

Since ijjPs defines a current, Pg defines a zero current situation or an equilibrium condition. 

The equilibrium solution Pg from (55) can now be used to calculate the steady state 
entropy production as given by Eq.(25). We thus have 

Explicit evaluation shows 

Sprod — As , (58) 

where As is given by Eq.(54). 
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The above result demonstrates a rather straightforward connection between entropy 
production-hke quantity of a chaotic system and the dynamics through dissipation constant 
7, parameters of the Hamihonian and correlation of fluctuations of the second derivatives 
of the Hamiltonian in the steady state. It is important to note that since one is concerned 
here with a few degree-of-freedom system with no explicit reservoir temperature does not 
appear in the expression for entropy production like term (58). The entropy production in 
a truly thermodynamic system and in the present case are therefore distinct. 



B. Numerical Verifications 



To verify the above theoretical analysis in terms of numerical experiments we now con- 
centrate on the following two points. First it is necessary to establish numerically that the 
dynamical system reaches a steady state, i.e., the probability density function P{{Xi},t) 
attains a steady state distribution Ps{{Xi}) in the long time limit. Second, the entropy pro- 
duction in the steady state calculated by the formula (58) needs to be verified numerically. 
To address the first issue we now proceed as follows ; 

The dissipative chaotic dynamics corresponding to the model Hamiltonian (39) is gov- 
erned by Eqs.(40) and (42). We choose the following values of the parameters |^ a = 0.5, 
b = 10, uj = 6.07 and 7 = 0.001. The coupling-cum-field strength e is varied from set to 
set. We fix the initial condition 2:1 (0) = —3.5, Z2{0) = which ensures strong global chaos 
P^ . To determine the steady state distribution function, say, Ps{Xi), where Xi = Aq (Eq. 
40) from the dynamical point of view we first define do as the separation of the two initially 
nearby trajectories and d{t) as the corresponding separation at time t. To express d{t) we 
write d{t) = [Ef (^i)^ + E?=Ar+i(^i)^]^- d{t) is determined by solving numerically Eqs.(40) 
and (42), simultaneously for the initial conditions of Zi and Z2 corresponding to Eq.(42). To 
follow the evolution of Xi numerically in time, i . e . , in going from j-th to (j+l)-th step 
of iteration, say, Xi has to be initialized as = ^d^. (The time evolution of the other 
components of X can be followed similarly). This initialization implies that at each step, 
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iteration starts with the same magnitude of do but the direction of do for step j+1 is that 
of d(t) for j-th step (considered in terms of the ratio -p). For a more pictorial illustration 



we refer to Fig. 1 of Ref. j-th term of iteration means t = jT {j = 1, 2, ■ ■ ■ oo), where 
T is the characteristic time which corresponds to the shortest ensemble averaged period of 
nonlinear dynamical system. 

Having obtained the time series in Xi (it may be noted that the series in Xj-s are also 

required for calculation of the largest Lyapunov exponent as defined by 

1 " d- 
X = lim In 

nT ^ do 

n ^ oo 
do ^ 

for the chaotic system) the stationary probability density function Ps{Xi) is computed as 
follows : The Xi axis ranging from -2 to +2 is divided into small intervals AXi of size 
0.025. The time series in Xi is computed over the time intervals of 1000-10000 times of 
the time period T. For each time interval AXi a counter is maintained and is initially set 
to zero before the simulation is started. The respective counter is incremented whenever 
Xi falls within the given interval. Finally the steady state probability distribution function 
Ps{Xi) is obtained by normalizing the counts. The result is shown in Fig.l for e = 10. 
Our numerical analysis shows that the distribution function attains stationarity at around 
t = lOOOT, beyond which no perceptible change in the distribution is obtained. 

We now turn to the second issue. In what follows we shall be concerned with the 
steady state entropy production (58) and its numerical verification. This quantity can be 
calculated in two different ways. First, it may be noted that the determination of Sprod 
(Eq.(57)) rests on two quantities defined in the tangent space; the steady state probability 
distribution function Ps{U) and the diffusion coefficient T>s in fZ-space. Once the procedure 
for calculation of the distribution function as illustrated above is known from the time 
series in Xi or X2 the evaluation of Ps{U) is quite straightforward since U is expressed 
as f/ = agXi + X2 according to (52). Here as is given by (54) with u"^ = 2b + c + C2 
and the average c and the integral over the correlation function, C2 are as defined in (51). 
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To calculate the correlation function {{Cit)((t — t))) and the average (C(^)) it is necessary 
to determine long time series in ({t) {C{t) = —12azf) by numerically solving the classical 
equation of motion (42) in phase space followed by averaging over the time series. For further 



details of numerical analysis of correlation functions we refer to the earlier work p!5HT7| . The 
diffusion coefficient can be numerically determined from the time series in U. Having 
known Ps{U) and T>s one can make use of the formula (57) to obtain the entropy production 
in the steady state. Sprod is thus numerically calculated. The second procedure of calculation 
of Sprod is the direct theoretical evaluation of As from the expression (54). Since the value 
of \s rests again on uj"^ and 7 and the dependence of uj"^ on the averages and the correlation 
functions are already known from the numerical analysis of phase space, can be calculated 
in the usual way. In Fig. 2 we compare the values of steady state entropy production thus 
obtained by the two different ways for several values of coupling-cum-field strength e. Here 
it should be noted that the curve connecting the squares (theoretically calculated entropy 
production in the steady state) corresponds to the negetive of entropy flux (Sjiux) since for 
the given example Sjiux in U space is — A^ (— / \sPs{U)dU, from Eq. (26)) for normalised 
probability distribution function Ps{U). Thus Fig. 2 is a numerical proof of Sprod = —Sfiux- 
The agreement is found to be quite satisfactory. We therefore conclude that at least for 
the model studied here and for the similar class of models the correspondence between the 
formulae of steady state entropy production and the numerical computation is fairly general. 



C. Entropy production in presence of weak forcing 

We now introduce an additional weak forcing in the dynamics. This is achieved by 
subjecting the dissipative chaotic system to a weak magnetic field (B) through a velocity 
{v) dependent force term ^^v x B where e and c; are the electric charge and the velocity of 
light. For simplicity we apply the constant field Bz which is perpendicular to gi-direction 
[see Eq.39]. In presence of this force field the motion of the particle will not be restricted to 
Bz and qi only. We shall have to consider the other direction q2 which is perpendicular to 
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both qi and Bz- 

To make the notation consistent with Eq.(8) we would now hke to correspond Xi,X2, 
and X4 to Aqi, Aq2, Api and Ap2, respectively. 

The relevant equations of motion are therefore as follows; 

X2 = X4 , 

Xs^2bXi-^Xs + C{t)X^ + h^Bz , (59) 

Q 

Xi = -h-fX^ - h^Bz . 

ci 

Here e/ci is the ratio of electric charge to the velocity of light used to put the equation in 
appropriate dimension. 

Then the non-equilibrium situation (due to additional forcing, h ^ 0) corresponding to 
Eq.(59) is governed by 

+ (60) 



where 



'Ds = V22 + T)2ias 



or more explicitly, 



2 



Proceeding as before we make use of the transformation of variables using 

U' = a'Xi + h'X2 + c'Xs + X4 (62) 
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so that Eq.(61) becomes 

dP d d'^P 

where 

P = P33c" + P3iaV . (64) 
and V^i in four dimension correspond to X>22 and 'D21 in two dimension, respectively 

and 

X'U' = -a'Xs - c'u^Xi + 7X3C' + h-fX4 - hB'^Xic' + hB'^Xs - hh'X^ . (65) 

Using (62) in (65) and comparing the coefficients of Xi-s we obtain 

6' = a' - c'7 + hB'^ - X'c' 

c' ^{h-f-X')/hB'^ (66) 

where A' is a solution of the cubic algebraic equation 

A'' + X'{h^^ + h^B'^^ - u^) - A''7(l + h)+ u^-fh = . (67) 

We now seek for a perturbative solution of the algebraic equation (67) which is given by 
(h as a small parameter) ; 

ft(A-g7-u-^-Ai7'-W;) 
° iX'l - 4Ai7 + A72 + _ ^2 

where Aq is the solution of (67) for h — 0; 

Ai = ^±4±^. (69) 

This is identical to A^ (Eq.54). Therefore by virtue of Eqs. (66-69) all the constants in (62), 
i.e., a', b', c' are now known. The stationary solution of (63) is now given by, 

Pl = N'e-^ , (70) 
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where N' is the normahzation constant. 

We are now in a position to calculate the steady state entropy flux ASfiux due to external 
forcing {h ^ 0) from Eq.(38) 

^Sfi^, = h'J dXdPV.F, + h'JdX ^1^^^) ■ (71) 

where the components of can be identified as 

Fii = Fi3 = B'^X4 
Fi2 = ^4 Fi4 = -B',X^ - 7X4 (72) 
and Vx-^i = —7 , 

and h5P = P'g — Pg denotes the deviation from the initial equilibrium state due to external 
forcing. For normalized probability functions and Pg the first integral in (71) vanishes. 
Since Pg is given by (55) with U as defined in (52), the expression for (71) reduces to 

= hB'^^ j (a,Xi + X3) dXP'g . (73) 

We now use the following transformations of variables 

u' = a'Xi + c'Xa + X4 (since h' = 0) 

v'^X^ 

w' = Xi 

and dX^dX^dX^^a'du'dv'dw' . (74) 
to calculate the integrals, 

V 



j X^XidXP'g 



2X'\c' 
V 



and Jx^X.dXP^^^ (75) 
which yield 
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'prod 



lux 




(76) 



For numerical verification of the above theoretical result (76) one can calculate entropy 
production (ASprod) numerically in the steady state in presence of weak forcing from Eq. 
(37) as in the previous subsection. The Eq. (37) for the present example reduces to the 
following form in the steady state, 



To calculate numerically ASprod ^^33,^^31 and Pg can be determined by directly using 
the procedure mentioned in the subsection VB and by simultaneously solving Eqs. (41) and 
(42). Similarly one can calculate Pg from Eqs. (42) and (59). Finally making use of all 
these quantities in Eq. (77) ASprod can be obtained. Thus the numerically evaluated ASprod 
should correspond to the results of Eq. (76) since our numerical verification in Fig. 2 shows 
good agreement between numerical and theoretical results, Dss^V^i being very close to X>22 
and V21, respectively, since h is very small. 

In the limit h and 7 small the above expression (76) can be simplified further. To this 
end we first note that 




(77) 



a^l ~ a; As ~ a; A' ~ Aq 



c = — and a = — is, . 
UJ 7 



(78) 



This reduces V further as follows; 



~ — (P33 - Vuu;) . 



(79) 



Thus we have 



v__ar_ 



(80) 
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Making use of (78-80), expression (76) can be approximated as 

^Sprod = hB'—- {—. + — 

uj^ 2 \a' a 



This expression is due to average of the work per unit time of the external force acting 
on the chaotic system. Note that the quadratic dependence on the magnetic field Bz in Eq. 
(81) is characteristic of an expression for entropy production in the steady state. Since the 
system is not thermostated this is independent of temperature. Although the leading order 
expression (81) is apparently free from diffusion coefficients, a close look into the more exact 
expression (76) reveals that their influence is quite significant in the higher order. 

VI. CONCLUSIONS 

Ever since the development of the theory of chaos, the dynamical variables in the strong 
chaotic regime have been interpreted as stochastic variables. One of the earliest well-known 



examples in this connection was set by demonstrating |TT[ the linear divergence of mean- 
square momentum in time in standard map, mimicking the Brownian motion. In the present 
paper we have tried to relate this chaotic diffusion to thermodynamic-like quantities by 
establishing a generalized Fokker-Planck equation pertaining to the tangent space. The 
explicit dependence of drift and diffusion terms on the dynamical characteristics of the 
phase space of the system is demonstrated. 

The main conclusions of our study are; 
(i) We analyze the nature of chaotic diffusion in terms of the properties of the phase space of 
chaotic systems. The drift and diffusion terms are dependent on the correlation of fluctua- 
tions of the linear stability matrix of the equation of motion. Since the later is the key point 
for understanding the stability of motion in a dynamical system, we emphasize that the 
thermodynamic-like quantities as discussed here have a deeper root in the intrinsic nature 
of motion of a few-degree-of-freedom system. 

25 



(ii) We identify the information entropy flux and production-like terms in the steady state 
which explicitly reveal their connection to dynamics through drift and diffusion term, in 
presence and absence of the external force field. 

(iii) The connection between the thermodynamically-inspired quantities and chaos are fairly 
general for the N-degree-of-freedom systems. 

The theory developed in this paper is based on the derivation of the Fokker-Planck 
equation for chaotic systems pertaining to the processes with correlation time which is short 
but finite (i. e. , for the systems with hard chaos). The suitable generahzation of the 
approach to more general cases, where one encounters long correlation time is worthwhile 
for further investigation in this direction. 
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APPENDIX A: THE DERIVATION OF THE FOKKER-PLANCK EQUATION 



We first note that the operator exp(— rVx ■ L°) provides the solution of the equation 
[Eq.®, a = 0] 

dfiX,t) 



dt 



-Vx-L7(X,t) 



(Al) 



/ signifies the "unperturbed" part of P which can be found exphcitly in terms of character- 
istic curves. The equation 

X = L°(X) (A2) 

determines for a fixed t a mapping from X{t = 0) to X{t), i. e., X — > with inverse 
(X^)-^ = X . The solution of Q is 

dx-' 



f{X,t) = fiX-\0) 



dX 



exp [-tVx-Fo]/(X,0). 



d(X) 



being a Jacobian determinant. The effect of exp(— tVx ■ L°) on /(X) is as 

dx-' 



exp(-tVx-L")/(X,0) = /(X"*,0) 



dX 



(A3) 



(A4) 



This simplification in Eq.(|T6D yields 

dP _ 
'dt ~ 



POO 


dX-^ 


/ dr 


Jo 


dX 



dX 



dX- 



P . 



(A5) 



{{Vx ■ L\X, t)Vx-r ■ Li(x-^ t - r))) 

Now to express the Jacobian, X^'^ and V^-^ in terms of Vx and X we solve Eq. (|A2| 
for short time (this is consistent with the assumption that the fiuctuations are rapid [0) 
Using Eqs.(§-^ we may rewrite "unperturbed" Eq. (|A2|) as 



d_ 
It 



^ X ^ 



and 



d 
dt 



X. 



2N 



( X ^ 



^2N 
( 



(A6) 



-7 
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+ 



G'Ar+i(X) 
G'2Ar(X) 



(A7) 



Here Gn+i{X) ■ ■ ■G2n{X) are the functions of {Xj} with i = 1- ■ - N only. This allows us 
to rewrite the solution of (|A6| ) and (^) as, 



+ 



and 



( y-r \ 



3 



y ^2iv y 



r 



y Gn{X) j 



(A8) 



(A9) 



G2n{X) 

Here the terms of O(r^) are neglected. Since the vector X~'^ is expressible as a function of 
X we write 



X'^ = G{X) 



(AlO) 



and the following simplification holds good; 



L\X-^, t - r).Vx- = L\G{X),t - r). Vx- 



d 



Y.Ll{G{X),t-T)^^ 



d 



Y: E Ll{G{X),t - r)9,k^ ; J, = 1 ■ ■ ■ 2X (All) 



j k 



where 



9jk 



dXr 



(A12) 



In view of Eqs.( |A8|) and ( |A9| ) we note: 



if j = k then gjk = 1, k = 1 ■ ■ ■ N 

= e-^\ k = N +1---2N 



if j ^ k then gjk oc —re 
or 
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Thus Qjk is a function of r only. 
Let 



From Eqs. 



Rj = Y.Ll{G{X),t-T)g,, 

k 

(g) and ( [MUD we write 

Ll{X-\t-T) = L]{G{X),t-T) =0 fori = 1 



..N 



(A13) 



(A14) 



So the conditions (|A13| ), ( |A14|) and ( [ASD imply that 



Rj{X, t-T) = Rj{Xi ■■■XN,t-r) for j = 1 • ■ ■ AT 
Rj{X, t-r) = Rj{Xi ■ ■ ■ X2N, t-T) ioi ] = N + 1---2N 



(A15) 



We next carry out the following simplifications of a^-term in Eq.(^). We make use of the 
relation (0) to obtain 

j j i * j j 



^L]{X,t)R, 



dXidXj 



P{X,t) 



where 



d 



(A16) 



i ^ 

The conditions ( A 14]) and (|A15| ) imply that 



(A17) 



R'j = for J = 1 ■ ■ ■ A^ 

R'j = Rj{Xi---XN,t-T) ^0; foTj = N + l---2N 



(A18) 



By ( [ATq ) one has 



R'.VxP{X,t) = Vx-R'P{X,t) 



(A19) 



Making use of Eqs.(0), (|ATTD , ( |ATB|) and (|ATU|) in Eq.(|A§) we obtain the Fokker-Planck 



equation ([T7|) 
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Figure Captions 



Fig.l. A plot of numerically calculated stationary distribution function Ps{Xi) as a function 
of Xi for the set of parameter values described in Sec.V. 

Fig. 2. The steady entropy production calculated numerically (circle) and theoretically using 
Eq.(58) (square) for different values of driven field strength e for the model described in 
Sec.V. 
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